Experimental observation of a strong mean flow induced by internal gravity waves 
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"^ , We report the experimental observation of a robust horizontal mean flow induced by 

^r^ \ internal gravity waves. A wave beam is forced at the lateral boundary of a tank flUed 

I— !■ with a linearly stratifled fluid initially at rest. After a transient regime, a strong 

p*^! jet appears in the wave beam, with horizontal recirculations outside the wave beam. 

^ ■ Using multiple scale analysis, we present a simple physical mechanism predicting the 
;— I \ 

. . growth rate of the mean flow and its initial spatial structure. We flnd good agreement 
c/3 . 

. y ' with experimental results. These results show that a mean flow with non-zero vertical 

^' vorticity can be generated by Reynolds stresses if the wave fulflls two conditions: i) 

, ^ ! the wave amplitude must vary along its propagation direction, which is the case in 

^r^ ' the presence of viscosity, ii) the wave amplitude must vary in the lateral direction, 

> ■ 

(^ ■ which is the case when the wave generator is localized in space. 
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Introduction. Stratified fluids support the existence of anisotropic dispersive waves, 
called internal gravity waves, which play a major role in astrophysical and geophysical 
fluid dynaniicsii^. Recent technical advances allowing for accurate visualization^i^ and well 
controlled wave generational^ in laboratory experiments have provided a renewal of interest 
in this fleld^. Previous laboratory experiments focused mostly on propagative wave beams 
in narrow tanks^'^'^i'^ or propagative vertical modesi^iii. Here we consider the case of a prop- 
agative wave beam in a wide tank, which remains largely unexplored, despite its physical 
importance. 

A central aspect of wave dynamics is the possible generation of a mean flow due to 
nonlinearities involving one or several wave beams. These phenomena have important con- 
sequences for geophysical flow modeling, since they imply backward energy transfers, or 
large scale transport properties induced by small scale motions, see e.g. Ref— >^. Among all 
waves, internal waves are very peculiar because of the speciflc and unusual nature of non- 
linearity. For instance, it has been reported^^ii^ that in some important cases, the leading 
nonlinear term unexpectedly cancels out if one has just one internal wave beam. 

King, Zhang and Swinney^^ have recently reported the generation of a mean flow by non- 
linearities in the presence of internal gravity waves. However, the structure of the observed 
mean flow has not been explained, and the underlying mechanism of generation has not yet 
been proposed. Lighthill^ noticed that either acoustic waves or internal gravity waves can 
generate steady streaming when their amplitude is attenuated by viscous effects. However, 
he did not performed explicit computation of the induced mean flow. Here we flU this gap 
by performing multiple scale analysis. We provide predictions for the spatial structure and 
temporal evolution of the mean flow which are supported by strong experimental evidences. 

There have been other theoretical studies of internal gravity waves-mean flow interactions, 
see e.g. Refs.— li^ii^, but none of them considered the case of propagative waves with a 
slowly varying amplitude in three dimensions. We will show that in our experiments, both 
viscous attenuation and lateral variations of the wave beam amplitude play a key role in the 
generation of the observed mean flow. 

The paper is organized as follows. We flrst present the experimental setup. Then we 
provide detailed observations of the wave fleld and of the mean flow. We flnally propose a 
mechanism to deduce the spatial structure of the mean flow and the temporal evolution of 
its amplitude from the measurements of the wave fleld. 



Experimental setup. We consider a 120 cm long, 80 cm wide and 42.5 cm deep wave tank, 
filled with 35 cm of salt water, see Fig. [U^a). The fluid is linearly stratified in density with 



a Brunt-Vaisala frequency A^ = \'—{g/p)dzP, where g is the local gravity, p the density 
of the fluid and z the vertical coordinate. An internal wave generator is placed on one 
side of the tank, see Refs.-"^ for its full presentation and characterization; it is made of a 
series of 18 rectangular plates stacked around a helical camshaft. The plates, which are 
14 cm wide in the y-direction, oscillate back and forth along the longitudinal horizontal 
coordinate x. The phase shifts between successive cams are chosen in order to form a 
sinusoidal profile at the surface of the wave generator. The rotation of the camshaft at a 
frequency u < N generates a moving boundary condition with an upward or downward phase 
velocity depending on the sign of the rotation of the helical camshaft. The displacement 
profile is Xo(t, z) = Xq sin(a;t — mz), with a vertical wavelength A^ = 27r/m = 3.8 cm and 
an amplitude Xq = 0.5 cm or Xq = 1 cm. 

In the experiments, the front face of the wave generator is located at a; = 0, centered 
at y = 0, z = 15.8 cm. The wave beam is L^h = 14 cm wide, 11.4 cm high, corresponding 
to three wavelengths. The wave generator is only forcing the x-component of the internal 
wave, while the z component is found to adjust according to the internal wave structure. 
The propagation angle 6 of the internal wave is varied by changing the rotation rate of the 
wave generator motor, while keeping the Brunt-Vaisala frequency constant between each 
experiment, namely A^ = 0.85 rad s"^. Importantly, this experimental set-up leads to a 
wave amplitude that depends on the frequency^. Moreover, the axis of the wave generator 
camshaft staying always vertical, the efficiency of the forcing depends significantly on the 
projection of the plate motion on the direction of propagation. The wave generator frequency 
being varied in the range u = 0.26 A^ to 0.50 A^, corresponding to an angle of propagation 9 
from 15° to 30°, the amplitude of the wave is measured experimentally. 

Velocity fields are obtained using a 2D particle image velocimetry (PIV) systent^. The 
fiow is seeded with 10 /ini tracer particles, and illuminated by a 532 nm 2W-continuous laser, 
shaped into a vertical or horizontal sheet. Respectively, a vertical 35x43 cni^ or horizontal 
33x43 cni^ field of view is acquired by a 8-bit 1024x1024 pixels camera. For each wave 
generator frequency, a set of 600 to 1600 images is recorded, at a frequency of 0.38 to 1.25 
Hz, representing 10 images per wave generator period. PIV computations are performed over 
successive images, on 21x21 pixels interrogation windows with 50% overlap. The spatial 



resolution is approximately 25 x 25 px/cm^. A snapshot of the particle flow is presented in 
Fig. [lb). 
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FIG. 1. (Color online) (a) Schematic representation of the experimental set-up. The wave generator 
is placed on one side of the tank, dehning the origin of the spatial coordinates. The excited 
plane internal wave has a frequency w, an upward phase velocity and propagates with an angle 
9 = sin~ {u)/N). (b) Top view of the particle flow in the horizontal plane z = 21.6 cm. 



Experimental results. Selective Fourier flltering of the measured horizontal velocity fleld 
allows to distinguish the "wave" part of the flow denoted by («', f ', w') from the "mean" part 
of the flow denoted by (m, tJ, w) . 

The wave part of the horizontal velocity fleld is obtained using a bandpass filter centered 
on the wave generator frequency w, with a width of 0.014 cj and is presented in Fig. [2](a) 
and (b). Figure [2](a) shows a vertical slice taken at the center of the wave generator (y = 0). 
We notice the transverse extension over three wavelengths, the amplitude decay in the x- 
direction. Figure [2](b) presents a horizontal slice located around the mid-depth of the wave 
generator showing the wave amplitude variations in the y-direction. The above experimental 
observations readily suggests that the wave part of the flow is monochromatic, propagating 
at an angle 9 with respect to the horizontal axis and with an amplitude varying slowly in 
space compared to its wavelength A: 



271" 

u' = U' cos \ ut ^ (2; sin 6 + z cos 6) 



(1) 



The knowledge of U' is sufficient to determine the other components f ', w' of the wave field 
oscillating at frequency a;, see Appendix. 
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FIG. 2. (Color online) (a) and (b): Experimental wave field, u', obtained by filtering the horizontal 
velocity field at frequency to. (c) and (d): Experimental mean flow, u, obtained by low-pass filtering 
the horizontal velocity field. The contours represent the amplitude of the wave horizontal velocity 
field, U'. The left panels (a) and (c) present the side view and the right panels (b) and (d) present 
the top view. All pictures were obtained for uj/N = 0.26 and a 1 cm eccentricity for the cams. 
The wave generator is represented in grey and the moving plates in black. The dashed line in (a) 
(respectively (b)) indicates the field of view of (b) (respectively (a)). 



Let us now consider the low-pass filtered flow, denoted as {u,v,w), which we call "mean 
flow". The width of the filter used to extract the mean flow is 0.25 a;. A strong jet going 
in the outward direction from the wave generator is observed in Fig. [2]^c) and (d). This 
structure is initially located close to the wave generator, then grows until it fills the whole 
plane. We observe that the jet is precisely produced inside the wave beam. Recirculations on 
the sides of the tank are visible in Fig. Wljo) and in Fig. El^d) through the blue patches (blue 
online). These recirculations clearly show that the horizontal mean flow is characterized by 
non-zero vertical vorticity defined as fi = d^v — dyU. Another striking observation is that the 



vertical velocity of the mean flow could not be distinguished from noise in the experiments. 
Consequently, the vertical component of the mean flow is negligible with respect to the 
horizontal components: w <^u. 

Interpretation. It has long been known that steady streaming associated with non-zero 
vertical vorticity can be generated by Reynolds stress due to an internal gravity (or acoustic) 
wave beam if the wave amplitude is attenuated by viscous effects^, but to our knowledge 
there exists no explicit computation of this phenomenon. Using multiple scale analysis, it 
is possible to go further and propose a prediction for the spatial structure and the temporal 
evolution of the mean-flow. 

Computations are presented in the Appendix, where we consider a slightly more idealized 
conflguration than our experimental setting in order to simplify the analytical treatment: 
the wave generator is assumed to be tilted in the direction of phase velocity and to be of 
infinite extension in this direction. We also assume that the velocity amplitude imposed by 
the wave generator varies smoothly from to [/ in the y-direction over a distance Ly^i,. The 
three independent adimensionalized numbers of the problem are i) the wave Froude number 
Fr^ = U / {NX) ii) the ratio X/Li, between wavelength and the attenuation length scale 
Lu = (A/ (27r)) N/u due to viscous damping in the direction of the wave propagationiii^, 
iii) the ratio \/ L^b between wavelength and the lateral extension of the wave generator. We 
consider a fluid in the Boussinesq approximatiort^ with the scaling FrA = e^, \/ L^ = e/Xy 
and \/ L^uib = ^/\i where \y and A^ are non-dimensional constants of order one and where e 
is the small parameter for the multiple scale analysis. In our experiments, these parameters 
are estimated as Frx ^ 0.2, X/L,^ ^ 0.2 and X/L^h ^ 0.2, so the assumptions underlying the 
multiple scale analysis are only marginally satisfled. However, this analysis will allow to get 
physical insight to the problem. The essential point is that nonlinear terms are negligible at 
lowest order, and that variations of the wave-amplitude in the y and x direction are of the 
same order of magnitude. 

A lengthy but straightforward computation presented in the Appendix shows that at 
order 5 the dynamics of the vertical component of the flltered vorticity is governed by 

dtU = (dy {w'd.u') - d^ {w'd.v')) + z/AH , (2) 

(corresponding to Eq. (12^ of the Appendix). Note that self-advection of the vertical 
vorticity by the horizontal mean flow remains negligible at this order. We see that nonlinear 
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terms act as a source of vertical vorticity, in the same manner as the Reynolds stress tensor in 
a turbulent flow acts as a source of turbulent transport. Using the ansatz ([1]) the nonlinear 
term of Eq. ([2]) can be exphcitly computed (see Appendix for more details) and Eq. ([2D 
reads: 

_ 8 r/'2 _ 

dt^ = "^ , + uAQ . (3) 

(2cos^)' ^ ^ 

A first consequence of Eq. ([3]) is that nonlinearities cannot be a source of vertical vortic- 
ity Q if the wave field is invariant in the y-direction. In our experiments, the wave generator 
occupies only one part of the tank width, so one might expect horizontal variations of the 
wave amplitude, which fulfills the necessary condition to observe generation of a horizontal 
mean flow associated with non-zero vertical vorticity. A second consequence is that if the 
wave field is symmetric in the y-direction then the source term in Eq. ([3]), resulting from 
a derivation of the wave field with respect to y, is antisymmetric in the y-direction. This 
is the case in our experiments, because the wave generator is symmetric with respect to y. 
A third consequence is that variations in the x directions are necessary to produce vertical 
vorticity. Such variation is due to viscous attenuationi>^ and occurs at a scale L^y introduced 
previously. 

Referencei^ reported that a single wave beam could not generate a mean flow in a two- 
dimensional (or axi-symmetric) configuration, with invariance in the y-direction. This is 
consistent with our finding that the generation of a mean flow requires variations of the 
wave-amplitude in the y direction. Refs.— "^ also showed that a mean flow can be generated 
by two interacting wave beams, even if there is no variation of the wave amplitude in the 
y direction. Again, there is no contradiction with our results, because the mean flow they 
describe is associated with zero vertical vorticity, due to their invariance in the y direction. 

Finally, let us describe the spatial structure of the nonlinear term of Eq. ([3]) in our 
experimental setup. The wave generator induces a symmetric wave field in the y-direction, 
maximum at the origin, satisfying {dxyU'"^) > for y > and {dxyU'"^) < for y < 0. 
According to Eq. ([3]), we conclude that the source term induces a dipolar vorticity structure 
associated with a horizontal jet going in the outward direction from the wave generator. 

Let us now estimate quantitatively this source term. In the experimental results presented 
here, the wave generator is not tilted with respect to the vertical and not infinite in the 
direction of phase propagation, but we expect that Eq. ([3]) remains valid if we consider 



the flow at mid-depth of the wave generator sufficiently far from the lateral boundary. 
Consistently with this remark, the horizontal slice shown in Fig. [2]^d) is centered on the 
middle plate of the wave generator (dashed line in Fig. [2]^ a)), For sufficiently short times, 
the viscous term of Eq. (^ can be neglected, and one expects that dt^ ^ {dxyU'^) / (2 cos 6) . 
The experimental determinations of these two terms are compared in Fig. [3] for t = 106 s. 
This time is much smaller than the viscous time (~ 1000 s) described in the last part of this 
paper. The magnitude of dtfl (Fig. [3t^a)) is of the order of the source term (Fig. [3](b)) and the 
two fields correspond well spatially. The inset in Figure [3]^a) shows a good spatial correlation 
between both fields. The slope that relates the two terms is 0.8 (±0.1), consistently with 
the theory, which predicts 1. The small discrepancy may be due the fact the wave generator 
is not tilted and not infinite in the direction of phase propagation. 
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FIG. 3. (a) and (b) present the top view in the z = 21.6 cm plane at t = 106 s of the two important 
quantities entering in Eq. ([3|): dt^ in the left panel and dxyU'^/ (2cos^) in the right panel. Both 
pictures were obtained for oj/N = 0.26 and a 1 cm eccentricity for the cams. The inset in the 
left panel presents the 2D-histogram of dtfl versus dxyW^ / {2 cos 6) to emphasize the correlation 
between both quantities. The dotted line corresponds to the slope 0.8 (±0.1) discussed in the text. 



In order to address the temporal evolution of the mean horizontal flow at mid-depth of 
the wave generator, it is convenient to integrate Eq. (^ over a domain defined as a half-plane 
of Fig.El^b) (domain V+ = [0,L^] x [0, Ly/2]). One obtains 

v+ ^ ^ v+ 



dxdy Ar2. 



(4) 



The integral in the left-hand side of this equation, which we call I{z, t), can be transformed 
as (if we denote C"^ the circulation path surrounding the domain P+) 

I(z,t) = \\ dxdy Q{x,y,z,t) = (p dl ■ u ^ / dx u, (5) 

^■^ Jc+ Jo 

since along C"^, the velocity is non zero only on the y=0 axis. In this form, I{z,t) can be 
identified as a measure of the jet strength. Equation (jlj) can then be written as 

Finally, defining Lq as the smallest characteristic scale, i.e. along y, of the vorticity fi, the 
viscous term may be approximated by —uI{z,t)/LQ. The jet strength is then solution of a 
first order differential equation, 

dJ{z,t) = S-^Iiz,t), (7) 

where S is the source term S = f/'^(0, 0, 2;)/ (2cos0) . Equation (jTj) shows that the jet 
strength I{z,t) should vary exponentially with time and one gets 

/(;,,t) = ^(^l_e-*Ag) . (8) 

We observe such an exponential growth in Fig. Hl^a), which represents the time evolution 
of the quantity I{z,t) for different values of S and z = 21.6 cm. A deviation from the 
exponential growth happens at large times when the mean flow reaches the side of the 
visualization window, in which case the integral I{z,t) saturates since the approximation of 
Eq. Ois not valid any more. The value of Lq is estimated by an exponential fit of the evolution 
of I{z, t) as a function of time. Remarkably, one gets the same characteristic scale Lq ~ 4 cm 
for all experiments, which corresponds to the jet width shown in Fig.[2t^d). We then estimate 
experimentally the source term in two different ways: the exponential fit (S'fit) and the 
measurement of the amplitude of the horizontal velocity field (^exp = {U'{0, 0, z)/ (2 cos 6*)) ). 
S'fit is plotted as a function of S^xp in Fig- IH^b). As expected, a linear relation between 
these two estimations is obtained. We find S'fit = (0.7 ± 0.1) Scxp in agreement with the 
spatial correlation observed in the inset of Fig. ^a). However our simple model predicts 
•S'fit = •S'cxp- We notice that the smaller the angle, the closer to the theoretical prediction is 
the experiment. The reason is the the wave generator is not tilted with the vertical axis in 



the direction of phase propagation, contrary to the assumption done for the multiple scale 
analysis, so our prediction become less accurate with increasing the angle 9. 

The fact that we obtain a fairly good prediction for the the long time evolution of the jet 
strength may seem surprising, since Eq. ([3]) holds only for sufficiently small time. The main 
reason is that the integration procedure over the half plane would cancel the additional 
horizontal self-advection term u// ■ VVt in Eq. (JSj), even if it may be locally important. 
Indeed, considering that the horizontal mean flow is non-divergent, the advection term can 
be written V ■ (^211//), and the surface integral of this term vanishes since Vt is zero on the 
Ox axis and the filtered horizontal velocity field is zero on the other edges of the domain. 
Consequently, we expect that Eq. ([8]) remains valid at time larger than those required to 
fulfill the hypothesis of the multiple scale analysis. 
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FIG. 4. (a) presents the evolution vs time of the quantity I(zo,t), with exponential fits, (b) shows 
the value S^t of the source term estimated from the exponential fits of panel (a) vs the experimental 
estimation Sgxp of the same source term. The dashed line corresponds to Sfit = (0.7 it 0.1) Sexp- 

Conclusion. We have reported experimental observations of a strong horizontal mean 
fiow with non-zero vertical vorticity when a propagative monochromatic wave is forced on 
the side of a tank filled with a linearly stratified fiuid. We stress here that there is no such 
mean fiow in absence of internal waves propagation, for example when the wave generator 
is excited at a frequency uo larger than the Brunt-Vaisala frequency A^, as in Ref.— . 

The key ingredient for the existence of this mean vertical vorticity with a dipolar struc- 
ture, associated with a strong horizontal jet fiowing outward of the wave generator is the 
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concomitant existence of variations of the wave amplitude in both horizontal directions. In 
the transverse direction (y), the variations are simply due to the fact that the wave genera- 
tor is localized in a segment smaller than the tank width. In the longitudinal direction, the 
variations of the wave amplitude are due to viscous attenuation. This shows the important 
role played by viscosity in the generation of the mean flow in our experiments. 

Interestingly, a similar phenomenon as the one described in this paper has been observed 
in the case of the reflexion of an internal wavebeam over topograph}^. Although viscosity 
is essential for the wave beam attenuation in the experiments, several other physical mech- 
anism may enhance its effect. For instance, parametric subharmonic instabilities followed 
by a direct energy cascade^i or wave breaking^^ may lead to a spatially varying wave beam 
amplitude. Rather than considering viscous effects, Ref.— models the effect of wave breaking 
by a linear damping term in the buoyancy equation and uses generalized Lagrangian-mean 
theory to describe the formation of a strong mean flow. This gives a complementary point 
of view to the present paper, which is focused on the description of a laminar experiment. 
We hope to address the effect of wave-breaking in future experimental work. 
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APPENDIX: MULTIPLE SCALE ANALYSIS 

We consider a Boussinesq fluid linearly stratified with buoyancy frequency A^ and viscos- 
ity z/, and a wave generator oscillating at frequency u = NsinO with a wavelength A. The 
wave generator is tilted at an angle 6 with the vertical axis. The amplitude of the horizontal 
velocity imposed by the wave generator is denoted as Ao{y)U. We adimensionalize length, 
time, velocity and buoyancy b by A/ (27r), 1/A^, U and NU, respectively. In terms of the 
velocity components [u, v, w), of the vertical vorticity Q = v^ — Uy and of the buoyancy b, 
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the dynamics reads 

nt + e^ {uh ■ Vh^ + (V^ ■UH)n + 9, {wd,v) - dy (wd^u)) = eX'^AQ, (9) 

Awu + Ahw = eX-^AAwt - e^ {dt (A^ (u • Vw) - d.Vn (u ■ Vuh)) + A^ (u ■ V6)) , (10) 

V/^ ■ u// = -u;^, (11) 

6i + e'^ {u-Vb)+w = 0, (12) 

where the indices t, x, y, z stand for partial derivatives, and where we have considered the 
scahng proposed in the text for the adimensional parameters (f// (AA^) = e^, vj (X'^N) = 
e/Xu). We denote uh = {u,v,0) the horizontal velocity field. Since the wave generator 
oscillates at frequency sin 6, we introduce the appropriate coordinates 



^ = X cos 6' — 2 sin 6' and r] = xsm9 + z cos 9 



(13) 



(see Fig. 5). The wave generator is located at ^ = and is of infinite extension in the 
//-direction. This amounts to impose the boundary condition 



wiC = 0, ?7, y, t) = Ao{y) sin (t sin 9 - r]) 



(14) 



The initial condition is such that there is no flow in the domain at t = {u = v = w = 
and fi = 0). We introduce the rescaled coordinates 



and 



^70 = V, 6 = eC, 6 = e^^ Vi = ey, U = e' with z = 0, 1, 



w = Wo (to, ti,t2, h, r], yi, ^i, ^2) + ewi (ti, ta, ^3, l/i, 6, 6) + o(e) 



(15) 



(16) 
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FIG. 5. The wave generator is assumed to be of infinite extension in the r/-direction, located at 
^ = 0. Arrows c,^ and Cg are phase and group velocities, respectively. 
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u = uo (to, ti, ta, ts, ?7, yi, ^i, 5) + ewi (ti, ta, Z/i, 6, 6) + o(e) (17) 

V = evi (to, ti, ta, ts, 77, yi, G, 5) + o(e) , (18) 

Q = e^Q2 {ts, 2/1,6,6) + e'fi4 (to, ti, ta, tg, r/, yi, 6, 6) + o(e') • (19) 

The strategy is the following: computations at order 0,1,2 allow to obtain the structure 
of the wavebeani by solving Eq. flTO|) . At these orders, nonlinear terms are negligible and 
there can consequently be no source of vertical vorticity. We then consider the dynamics of 
the vertical vorticity given by Eq. (jH]) at order 3,4,5 to find the structure of the vortical 
flow induced by the wavebeam through Reynolds stresses, considering a temporal filtering 
over one wave-period. Our scaling is such that generation of a mean-flow through Reynolds 
stresses due to the wave attenuation appears at order 5, and requires the knowledge of the 
wave structure at order 2. Note that using the notations introduced in the text, the mean 
vortical flow is given hj Q = Q2 = ^2 while the wave field is {u',v',w') = {uQ,evi,Wo). 
Indeed, we can show a posteriori that ^24 = at order 5. 

At zeroth order, we recover from Eq. (fTO!) the inviscid dispersion relation, while the boundary 
condition imposes the wavelength of the wave: 

Wo = 1^06-''"'+'*° ^'"^ + W*e'''°-''° ^'"^ + U, Urjo = . (20) 

Equation (ITT]) gives the corresponding velocity in the x-direction UQrjo = WQrjg / ta.n 6 . Con- 
sidering order 1 terms in Eq. ( ITOl) . we obtain that for ti > ^i/cos6', 

iyo = Ae-«^/(2"''^°^''), with A(t2,t3,yi,6 = 0) = Ao(yi), (21) 

Using Vlx^;, = Uoy^ leads to 



icosO 



Vl 



-5i/(2Ay cose) / yi* irjo-itosin 



.2'^ V-^Oj/i 



(^*^^g.,o-itosine _ J^^^^^-^,o+^tosine^^ + f^^ f^^^ ^ g . (22) 



(^smi ^ 

Considering the order 2 terms in Eq. (1101) . one can show that for ta > 6/ cos 6', Wq does not 
depend on ta. At order 3, all nonlinear terms appearing in Eq. ([9]) for the dynamics of Vt 
vanish. At order 4, we obtain 

^4*0 + ("^xo {wod^^vi) - dy^ {wodzaUo)) = . (23) 

Using Eqs. (I2m22p and moj^q — u!0'no/ '^^'^d, one can show that d^^^ (wodz^^Vi) — dy^ (wodzoUo) 
oscillates at frequency 2 sin 9 and that the term independent of tg vanishes: the nonlinear 
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terms do not induce a steady streaming at this order. 

We introduce the filtering over one wave period fi = (sin6'/27r) L'^ dt fi. This definition 

is shghtly different from the filtering procedure in our experiment, but we keep the same 

notation. Filtering Eq. ([9]) at order 5, and using the results obtained at order 0-4, we find 

that 

^2*3 = [dyi (Wodz^Uo) - 9^1 (WodzoVl)] + X~^ (^26a + ^-iyiyi) , (24) 

with Q2 = f22-Using Eq. (I2m22p as well as Wor^o = Worj^/tanO, the Reynolds stresses can be 
explicitely computed. For ti > .^i/cos6', ^2 > 'C2/COS6', One finds that 

O _ ~^ 9y^\\^o\ ) -Ci/(X.cos9) I 1 /Q , Q ^ (0K\ 

^hts - . 2 ^^ ^ + IT '^^^2651 + ^hyiyi) {^^) 

(Sm &) Aj, A,y 

Experimentally, we measure the wave amplitude of the velocity in the x-direction U' = 

2\Wo\/tand. Using Eq. ( l2Tll22|) and Wor^o = w^Orjo/^^'^^' -^o can be expressed in term of U'. 

Coming back to dimensional units we obtain the following equation for the temporally filtered 

vorticity 

n, = ^^1^ + uAU . (26) 

(2cos^)' ^ ^ 
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